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Line formation in convective stellar atmospheres 

I. Granulation corrections for solar photospheric abundances 
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Abstract. In an effort to estimate the largely unknown effects of photospheric temperature fluctuations on spectroscopic abun- 
dance determinations, we have studied the problem of LTE line formation in the inhomogeneous solar photosphere based on 
detailed 2-dimensional radiation hydrodynamics simulations of the convective surface layers of the Sun. By means of a strictly 
differential 1D/2D comparison of the emergent equivalent widths, we have derived "granulation abundance corrections" for 
individual lines, which have to be applied to standard abundance determinations based on homogeneous ID model atmospheres 
in order to correct for the influence of the photospheric temperature fluctuations. In general, we find a line strengthening in the 
presence of temperature inhomogeneities as a consequence of the non-linear temperature dependence of the line opacity. The 
resulting corrections are negligible for lines with an excitation potential around Ei = 5 eV, regardless of element and ionization 
stage. Moderate granulation effects (A gran ~ —0.1 dex) are obtained for weak, high-excitation lines (Ei J> 10 eV) of Cl, Nl, 
Ol as well as Mgll and Sill. The largest corrections are found for ground state lines (Ei = eV) of neutral atoms with an 
ionization potential between 6 and 8 eV like Mg I, Ca I, Ti I, Fe I, amounting to A gran ~ — 0.3 dex in the case of Ti I. For many 
lines of practical relevance, the magnitude of the abundance correction may be estimated from interpolation in the tables and 
graphs provided with this paper. The application of abundance corrections may often be an acceptable alternative to a detailed 
fitting of individual line profiles based on hydrodynamical simulations. The present study should be helpful in providing upper 
bounds for possible errors of spectroscopic abundance analyses, and for identifying spectral lines which are least sensitive to 
the influence of photospheric temperature inhomogeneities. 
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1. Introduction 

Theoretical model atmospheres serve as the fundamental tool 
for the quantitative analysis of stellar spectra. The basic param- 
eters of a stellar atmosphere - effective temperature, surface 
gravity, and chemical composition - can be inferred from the 
comparison between observed and calculated spectra. Much ef- 
fort has been spent to improve the theoretical models by incor- 
porating as much realistic physics as possible. Specifically, in 
late-type stars one is confronted with the problem of a detailed 
treatment of both radiative and convective energy transport. A 
high degree of sophistication has been achieved in modeling 
radiative transfer: standard stellar atmosphere codes can now 
handle the influence of millions of spectral lines. 

In contrast, the treatment of convection is still rather crude. 
Usually, stellar atmospheres are described as 1 -dimensional hy- 
drostatic configurations where the convective energy transport 
is calculated from the so-called mixing-length theory (MLT, 
Vitense 1953|; Boh m-Vitense 1958) or variants thereof (e.g. 
Canuto et al.1996). However, modern spectroscopic observa- 



tions have reached a level of precision that requires improved 
model atmospheres for a reliable interpretation, including a 
proper treatment of hydrodynamical phenomena. 

From the point of view of standard model atmospheres, 
convection affects the temperature structure in a twofold way: 
it influences the mean vertical stratification and introduces hor- 
izontal inhomogeneities. MLT is designed to model only the 
mean structure and is not suitable to construct realistic multi- 
component model atmospheres. Although realistic hydrody- 
namical model atmospheres do exist for the Sun and a few 



solar-type stars due to the pioneering work by Nordlund (1982) 
Nordlund & Dravins (119901), Stein & Nordlund (|1989l 1998) 
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and Ludwig et al. (1999), to name just a few important contri- 
butions, such complex models are not usually applied for the 
analysis of stellar spectra. 

Observationally, the solar granulation is the visible imprint 
of thermal convection extending far into the photospheric lay- 
ers where spectral lines are formed. The effect of the associated 
photospheric temperature inhomogeneities on the line forma- 
tion process and the consequences for spectroscopic abundance 
determinations are not yet well understood, but it is generally 
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believed that the errors introduced by representing a dynamic, 
inhomogeneous atmosphere by a static, plane-parallel model 
are small. However, reliable quantitative estimates of this ef- 
fect are difficult to find. Earlier studies suffer from uncertain- 
ties in the underlying empirical two-component models, which 
are rather ad hoc and lack a firm physical foundation (e.g. 
Hermsen 1982). Later investigations based on relatively crude 
numerical convection models (Atroshchenko & Gadun 1994; 



Gadun & Pavlenko 1997) led to questionable and partly am- 
biguous conclusions. Recent determinations of the solar pho- 



tospheric Fe and Si abundance by Asplund et al. (2000b) and 



Asplund (2000a), respectively, rely on state-of-the-art 3D nu- 



merical convection simulations and take all hydrodynamical ef- 
fects fully into account. However, comparison of these results 



with those based on the empirical ID Holweger-Miiller (1974) 



solar atmosphere does not allow to separate the influence of in- 
homogeneities from the effect of the mean stratification, and 
to study the 1D/3D abundance corrections as a function of the 
line parameters in a systematic way. 

In this work we investigate the spectroscopic effects of hor- 
izontal temperature inhomogeneities in the solar atmosphere 
based on detailed 2D hydr odyna mical models o f solar surface 
conve ction (Freytag et al. 1996| ; Ludwig et al. 1999 , Steffen 
200C ). We do not address the question whether the mean tem- 
perature structure supplied by MLT is an appropriate represen- 
tation of the mean stratification obtained from hydrodynamical 
simulations (the problem of the 'right' choice of the mixing- 
length parameter a for recovering the correct mean tempera- 
ture stratification was investigated by Steffen et al. (1995 ) in the 
context of white dwarf atmospheres). Rather, our main question 
is: how large are the systematic errors of standard spectroscopic 
abundance determinations due to the fact that one replaces the 
'real' line profile, which actually is the result of averaging the 
spatially resolved line intensities over the granulation pattern, 
by the line profile computed from a single plane-parallel tem- 
perature stratification? 

In Sect. |2|, we introduce a simple analytical argument to il- 
lustrate the basic effect of line strengthening in the presence of 
temperature fluctuations. Then we summarize the main char- 
acteristics of the numerical simulations of solar surface con- 
vection in Sect. ||, outline our method to derive 'granulation 
abundance corrections' from a differential 1D/2D comparison 
of computed line equivalent widths in Sect, [ij before finally 
presenting the results for the Sun in Sect. ^. A summary of our 
findings is given in Sect. ^. 

In order to give a physical explanation for the numerical 
results obtained in this work, we have carried out a detailed 
investigation the process of LTE line formation in an inhomo- 
geneous stellar atmosphere based on the analysis of the trans- 
fer equation and the evaluation of line depression contribution 
functions. This study is presented in the second paper of this 
series (Paper II). 

2. A simple analytical illustration 

The strength of a spectral line formed in a plane-parallel (ID) 
atmosphere will differ from the strength of the same line 
formed in an inhomogeneous atmosphere with the same chemi- 




Fig. 1. Spatial variation of W for a sinusoidal temperature vari- 
ation according to Eq.(^), for a = 0.05 and Ei = 10 eV. 
Averaging over one spatial wavelength gives (W) /Wo = 1.78 
(dotted horizontal line). 



cal composition and mean temperature stratification because, in 
general, there is a non-linear relationship between line strength 
and temperature. It may be expected that lines with a very high 
excitation potential are most affected because (i) they orig- 
inate from deep layers where convection is well developed, 
and (ii) they are highly sensitive to temperature fluctuations. 
Likewise, the molecular equilibrium is known to be highly sen- 
sitive to temperature fluctuations, so molecular lines might also 
be strongly affected if the temperature inhomogeneities reach 
the higher photosphere where these lines are formed. 

As an illustration of the expected effect, consider the fol- 
lowing simple analytical argument. Assume that the equivalent 
width of a weak line W depends exponentially on the inverse 
temperature 8 — 5040/T, as suggested by the form of the 
Saha-Boltzmann equations, [] 



W(A6) = Wq exp {-E % A9 In 10} . 



(1) 



where Ei is the excitation potential in eV and AO = 9 — 6q- 
Assume further that 9 is normally distributed about the mean 
9 with standard deviation er, 



1 [ A6 2 



(2) 



Then it can be shown by analytical integration that 

/+oo 
g(A9)W(A9) dA9 
-OO 



Wo exp <{ |£ 2 cx 2 (lnl0) 2 



(3) 



i.e. the mean value (W) of the equivalent width formed in an 
inhomogeneous medium is always greater than the equivalent 
width for the unperturbed medium, Wo- As long as EiU is small 
enough, the effect is negligible, but this is not true in general. 
Assuming a = 0.05, corresponding to 6% temperature fluctu- 
ation around Tq = 6000 K, and Ei = 10 eV, the effect is very 



1 Note that Eq.([j]) is a reasonable approximation only for weak lines 
of 'majority species' with high excitation potential (cf. Gray 1992). 
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substantial: (W) » 2 Wq, i.e. the elemental abundance derived 
from such lines would be overestimated by a factor of 2! Note 
that, for example, several important neutral oxygen lines have 
Ei between 9 and 1 1 eV. 

A similar result is obtained when assuming a totally differ- 
ent distribution of Ad, 



/(A0) 



1 



7rV2a 2 - A6» 2 ' 



(4) 



corresponding to sinusoidal temperature fluctuations in one 
horizontal direction, 



A0 = 0"\/2sin(27nE), 



(5) 



where x is the spatial coordinate. According to Eq.(|lJ), the fluc- 
tuation of the equivalent width is then 

W{x) = Woexp|-EilnlO(rV2sin(2TO)| . (6) 

Adopting again a = 0.05 and Ei = 10 eV, we have plotted 
W(x) in Fig. [jj Clearly, the fluctuation of W is highly asym- 
metric with respect to Wq, and numerical averaging over one 
spatial wavelength gives (W) ~ 1.8 Wo- 

Looking at Table [j], it turns out that the simple model 
severely overestimates the abundance correction in the case of 
the oxygen lines mentioned as an example above. This is be- 
cause, in reality, the line formation process is much more com- 
plicated, since the line strength depends not only on the line 
opacity but also on the continuum opacity and on the local tem- 
perature gradient of the atmosphere (see Paper II). Moreover, 
spectral lines usually form over a substantial depth range, and 
the sign of A9 may depend on depth, leading to partial cancel- 
lation of the local effects, especially for stronger lines. Finally, 
saturation can play a significant role as well. Hence, realistic 
abundance corrections can only be derived from detailed line 
formation calculations based on multi-dimensional hydrody- 
namical model atmospheres. 

3. Time-dependent Radiation Hydrodynamics 

In order to adequately address the complex problem outlined 
above, we have investigated the effects of temperature inho- 
mogeneities on the formation of a variety of spectral lines by 
means of a realistic 2-D radiation hydrodynamics (RHD) sim- 
ulation of solar surface convection. In the following, we briefly 
summarize the foundations of the hydrodynamical simulation 
and the main results concerning the solar photosphere. 

3. 1. 2-dimensional numerical simulation 

Stellar surface convection is governed by the conservation 
equations of hydrodynamics, coupled with the equations of 
radiative energy transfer. Our hydrodynamical models result 
from the numerical integration of this set of partial differen- 
tial equations. This approach constitutes an increasingly pow- 
erful tool to study in detail the time-dependent hydrodynamical 
properties of the solar granulation. 

Just like 'classical' stellar atmospheres, the hydrodynami- 
cal models are characterized by effective temperature, T e g , sur- 
face gravity, g, and chemical composition of the stellar matter. 



But in contrast to the mixing-length models, they account for 
'overshoot' and no longer have any free parameter to adjust 
the efficiency of the convective energy transport. Based on first 
principles, RHD simulations provide physically consistent ab 
initio models of stellar convection which can serve to address 
a variety of questions, including the problem of line formation 
in inhomogeneous stellar atmospheres. 

The models used for the present investigation comprise a 
small section near the solar surface, extending over 7 pres- 
sure scale heights in the vertical direction, including the photo- 
sphere, the thermal boundary layer near optical depth tr oss = 
1, and parts of the subphotospheric layers. Only the upper- 
most layers of the deep solar convection zone can be included 
in the model, requiring an open lower boundary. The simula- 
tions are designed to resolve the solar granulation. The effect 
of the smaller scales, which cannot be resolved numerically, 
is modeled by means of a subgrid scale viscosity (so-called 
Large Eddy Simulation approach). Spatial scales larger than 
the computational box are unaccounted for. We employ a re- 
alistic equation of state, including ionization of H, Hel, Hell 
and H2 molecule formation. In order to avoid problematic sim- 
plifications like the diffusion or Eddington approximation, we 
solve the non-local radiative transfer problem along a large 
number (iV x / sub N e N$ = 210 • 4 • 2 • 8 « 10 000) of rays in 
5 opacity bins, with realistic opacities accounti ng ad equately 
for th e influence of spectral lines (see Nordlund 1982 , Ludwig 



1992[). Fo r further details see L udwig et al. (|1994j ), Freytag et 



al. ( |1996| ), Ludwig et al. ( |1999| ), Steffen (gOOOj) 



3.2. Main results 

The basic features of the numerical convection model are illus- 
trated in Fig. ||, showing a representative snapshot from a well- 
relaxed simulation of the solar granulation. In the convectively 
unstable, subphotospheric layers, fast (v c s ) narrow down- 
drafts (also called plumes or jets) stand out as the most promi- 
nent feature. They are embedded in broad ascending regions 
(the granules) where velocities are significantly lower. We find 
enormous temperature differences between the cool downflows 
and the hot upflows (AT max « 5 000 K at z m -100 km). 
Near optical depth t = 1, efficient radiative surface cooling 
produces a very thin thermal boundary layer over the ascend- 
ing parts of the flow, where the temperature drops sharply with 
height (AT/ Az w 80 K/km, see upper panel of Fig. |J) and 
the gas density exhibits a local inversion. Note however, that 
the temperature contrast is much reduced on the optical depth 
scale (see bottom panel of Fig. |]). 

Although convectively stable according to the 
Schwarzschild criterion, the photosphere (r < 1) is by 
no means static. Convective flows overshooting into the 
stable layers from below are decelerated here and deflected 
sideways. This can sometimes result in transonic horizontal 
streams which lead to the formation of shocks in the vicinity 
of strong downdrafts (see vertical fronts in Fig. ^| between 
3 000 < x < 3 500 km, z > 100 km). 

Oscillations excited by the stochastic convective mo- 
tions contribute to the photospheric velocity field as well. 
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Solar Granulation: Model L71 D09, t= 41 940 s (TEMPERATURE) m<p) iog M (t) 




1000 2000 3000 4000 5000 

x [km] 



Fig. 2. Snapshot from a 2-dimensional numerical simulation of solar surface convection after 41 940 s of simulated time. This 
model was computed on a Cartesian grid with 210x106 mesh points (tick marks along upper and right side), with periodic lateral 
boundary conditions (L = 5 250 km). "Open" boundary conditions at the bottom and top of the computational domain are 
designed to minimize artificial distortions of the flow. The velocity field is represented by pseudo streamlines, indicating the 
displacement of a test particle over 20 sec (maximum velocity is 10.1 km/s at this moment); the temperature structure is outlined 
by temperature contours in steps of 500 K. Geometrical height z = (scale at left) corresponds to tr oss w 1; scales at right refer 
to the horizontally averaged gas pressure [dyn cm~ 2 ] and Rosseland optical depth tr oss . 



Interestingly, the oscillation periods lie in the range 150 to 
500 s for the solar simulation, in close agreement with the 
observed 5 minute oscillations. Controlled by the balance be- 
tween dynamical cooling due to adiabatic expansion and radia- 
tive heating in the spectral lines, the mean photospheric temper- 
ature stratification is slightly cooler than in radiative equilib- 
rium, and the resulting temperature structure is not at all plane- 
parallel. 

According to Eq.(|3|), the line strengthening depends on 
the key quantity a = 9o <5T rms /Xo, where ST rms = 
yj ((T — To) 2 ). We have evaluated <5T rms as a function of op- 
tical depth for our 2D numerical simulation of solar convec- 
tion (Fig. Q). For the problem under investigation, tempera- 
ture differences have to be taken over levels of constant opti- 
cal depth, rather than at constant geometrical depth. According 
to the simulations, typical values in the lower photosphere 
are in the range 0.04 < ST lms < 0.06. The minimum near 
logTRoss = —0.5 is related to the 'inversion' of the tempera- 
ture fluctuations at this height. 

It is worth mentioning that we have recently carried out 
fully 3-dimensional simulations with a new code developed by 
B. Freytag and M. Steffen (Freytag; Steffen et al., in prepara- 
tion). We find that the depth dependence of ST Tms (r) is qual- 
itatively similar in 2D and 3D. However, the amplitude of the 
temperature fluctuations is systematically smaller at all heights 
in the 3D simulation. The difference is most pronounced in the 
higher photosphere (see Fig. Q). 

3.3. Accuracy of 2D and 3D convection models 

The straightforward application of hydrodynamics and radia- 
tive transfer simulations is certainly the most satisfactory and 
consistent method to treat the line formation in a convec- 
tive stellar atmosphere, although the outcome is difficult to 
check by independent analysis. The hydrodynamical models 



presented here must be regarded as somewhat preliminary for 
several reasons. First of all, the restriction of the flow geome- 
try to 2 dimensions must lead to deviations from the true 3D 
atmospheric structure. As a consequence, the convective line 
shift and asymmetry predicted by the 2D granulation models 
is quantitativel y differ ent from that obtained with 3D models 
(Asplund et al. 2000a I. A further problem is the relatively poor 
spatial resolution of the simulations (especially in the vertical 
direction), which is barely sufficient to resolve the thin thermal 
boundary layer near t = 1. Certainly, the Reynolds number 
of the simulations is many orders of magnitude below realistic 
values. Finally, the influence of the spectral lines on the ra- 
diative energy balance can be treated only in an approximate 
way. Typically, a few opacity bins are used in the simulations, 
as compared to thousands of frequency points in standard ID 
models. 

The latter problems apply to state-of-the-art 3D convec- 
tion simulations as well. We note that a perfect line profile fit- 
ting does not necessarily imply a perfect model atmosphere. 
Another important diagnostic is the center-to-limb variation of 
the continuum intensity. Even the best 3D solar granulation 
models still fail to reproduce the observed continuum proper- 



ties with high accuracy (Asplund et al. 1999), indicating that 



"the temperature structure close to the continuum forming lay- 
ers could be somewhat too steep". 

In view of the above mentioned limitations, we prefer 
to avoid a direct application of our convection simulations. 
Instead, we favor the somewhat less powerful but more reliable 
differential approach described in the following section. 

4. Granulation abundance corrections 

4.1. The problem 

The problem we want to address is the following. Assume 
you have the 'best possible' ID representation of the solar at- 
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Stellar granulation model L71 D09_625x, Time= 41 940 s 




-1000 -500 

Height z [km] 

Stellar granulation model L71 D09_625x, Time= 41 940 s 



500 



T(i) snap average 
T (t) snap average 




Fig. 3. Temperature as a function of geometrical depth z (top) 
and as a function of individual optical depth tr oss (bottom) 
for the snapshot shown in Fig. ||. Each line corresponds to a 
different x position. Note that the (T) and the (T 4 ) 1 / 4 average 
are hardly distinguishable in this representation. 



mosphere, correctly reproducing the continuum intensity dis- 
tribution I\ (A) and its center-to-limb variation I\{fi). The 
microturbulence parameter £ mic is adjusted such as to mini- 
mize the overall dependence of the derived abund ances on line 
strength. The empirical Holweger-Muller ( 1974 ) atmosphere 
comes close to this ideal mean temperature structure. In the 
presence of temperature fluctuations, however, non-linearities 
in the line formation process will inevitably lead to systematic 
errors of abundance determinations based on individual spec- 
tral lines. The present work aims at quantifying this kind of 
error. 

We have adopted a strictly differential procedure to quan- 
tify these errors. The basic idea is to identify our 2D hydrody- 
namical convection simulations with the real inhomogeneous 
atmosphere, and the mean structure of the simulations with the 
'best possible' ID representation. Note that it is not impor- 
tant for the hydrodynamical convection models to reproduce 
the mean temperature structure of the real solar atmosphere 
with great precision, because the related errors cancel to a first 
approximation in our differential comparison. Similarly, the 
details of the hydrodynamical velocity field are of secondary 
importance in the present context, since such 1D/2D differ- 




4 2 -2-4 

Optical Depth Log 10 (T Ross ) 

Fig. 4. Time-averaged amplitude of rms temperature fluctua- 
tions as a function of Rosseland optical depth, derived from the 
numerical simulation of solar granulation described in Sect. |3| 
(solid). r5T rms was obtained by averaging over surfaces of con- 
stant Rosseland optical depth for each snapshot and subse- 
quently taking a long-term time average. For comparison, we 
show corresponding results from a 3-dimensional numerical 
simulations (dashed) recently carried out with a new code by 
B. Freytag and M. Steffen. 

ences are largely eliminated by a proper choice of the micro- 
turbulence parameter. In contrast, there is no free parameter 
to account for the presence photospheric temperature inhomo- 
geneities. This is why we focus on the abundance corrections 
related to temperature fluctuations <5T rms only. The most impor- 
tant quantity to be provided by the convection model is there- 
fore the magnitude and height dependence of 8T Yms . 

In the following, we outline our method to derive LTE 
"granulation abundance corrections", A glan , which must be 
added to the logarithmic elemental abundance derived from 
standard ID model atmosphere analysis in order to compen- 
sate for the effects of temperature inhomogeneities. Note that 
our procedure ensures that the corrections A gran — > for van- 
ishing horizontal temperature and pressure fluctuations. 



4.2. The procedure 

First we have selected S = 11 representative snapshots from 
our time-dependent convection simulation (see Sect. |J). These 
snapshots (one of them shown in Figs. || and |]) all have the 
correct effective temperature (the horizontally averaged radia- 
tive energy flux lies within ±1% of the nominal solar flux), but 
differ in the number and size distribution of granules and in the 
continuum intensity contrast (SI lms varies between 18.0 and 
26.5% at A 5000 A). For each of the 2D snapshots (s = 1 . . 11) 
we then constructed a corresponding ID model by averaging 
temperature T and gas pressure P over the corrugated surfaces 
of constant Rosseland optical depth: 



TmdItross) = (^ s 4 2d( t Ros S ,2;)) 1/4 , 

-Ps,lD(TRoss) = (-Ps,2D(TRoss, x)) , 



(7) 
(8) 
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where (.) denotes averaging over the horizontal coordinate x. 
Note that the ID models constructed in this way have very 
nearly the same radiative flux as the underlying 2D models. 

Next we compute, for given elemental abundance and 
atomic line parameters, the equivalent width of any given spec- 
tral line from (i) the set of 2D snapshots by horizontal averag- 
ing of the spatially resolved synthetic line profile 



W S .2V = (Ic, B (a:) W s>2 d(x)) / (ic,s(z)> 



(9) 



where I c {x) is the local continuum intensity, and (ii) from the 
corresponding ID models, resulting in W Si id- 

For the line formation calculations, we employed the Kiel 
spectrum synthesis code LINFOR under the assumption of 
LTE. For the ID atmospheres, we have always adopted a depth- 
independent micro turbulence velocity of £ m ; c = 1 km/s. For 
snapshots from the convection simulations we have computed 
synthetic spectra for two different assumptions concerning the 
non-thermal velocity field: in case (A) we used the depth- 
dependent hydrodynamical velocity field z) obtained from 
the simulations, in case (B) we replaced the hydrodynamical 
velocity field by a depth-independent microturbulence velocity 
Cmk = 1 km/s. 

As the final step, we calculate the mean equivalent widths 

resulting from the set of ID and 2D line profiles, Wid = 

i s i s 

s E s =i Ws,id and W 2 t> = $ Z)a=i W *,2D, respectively. For 
truly weak lines, the (logarithmic) "granulation abundance cor- 
rection" would simply be 



v wcak lines 
gran 



= logWi D -logW2D. 



(10) 



In general, however, the equivalent width is not propor- 
tional to the elemental abundance due to saturation. We ac- 
count for these curve-of-growth effects in the following way. 
Knowing that the 'true' abundance, A 2 d produces the equiva- 
lent width W 2 d with the 2D inhomogeneous model atmosphere 
(and Win with the ID mean structure), we iteratively deter- 
mine the elemental abundance Aid which produces the equiv- 
alent width W2D with the ID mean stratification. The "granu- 
lation abundance correction" is then given by 



log A 2D - log Ai D . 



(11) 



For simplicity, we u sed th e standard HM model atmosphere 
(Holweger & Miiller |1974| ; £ mlc = 1 km/s) for this final step 
of the procedure (Win ^ A 2 d, Wzd — * Aid). This is justi- 
fied because all we need is a realistic curve-of-growth, W(N), 
which is very similar for the averaged simulation snapshots and 
the HM atmosphere. 

Case (A) and (B) give essentially the same corrections for 
weak lines, but the results may differ substantially for stronger 
lines, where W2D.A < W2D.B (implying that the hydrody- 
namical velocity field provides an effective 'microturbulence' 
£mic(RHD) < 1 km/s). For lines of arbitrary strength, we 
strongly favor the granulation corrections derived from case 
(B), because this strictly differential 1D/2D comparison is 
based on the same non-thermal velocity field, and thus cleanly 
separates the effect of the temperature inhomogeneities from 
that of the small-scale velocity field. In the following, A gran 



therefore denotes the granulation abundance corrections de- 
rived from case (B). 

For practical reasons, we have used only a small number of 
snapshots (S — 11) for our present investigation. The statistical 
uncertainty of the derived abundance corrections is therefore 
appreciable. In order to quantify this uncertainty, we have also 
tabulated the standard deviation of the corrections A„ r!1T1 , 



\ 



1 6 
S- 1 5 



(12) 



derived from the 1 1 individual snapshots. The quantity a can 
serve to judge the significance of the abundance corrections for 
each of the lines listed in the tables below. 



5. Results for the Sun 

We have applied the above procedure to a representative sam- 
ple of lines, covering different elements, ionization stages, ion- 
ization energies, excitation energies, line strengths and wave- 
lengths. In this differential study, we have adopted the same 
representative van der Waals broadening constant for all lines, 
logC 6 = -30.3867 (definition accor ding t o Mihalas ( |1978[ ), 



which differs from that by Unsold ( 1955| ) by a factor 2n), 
and have neglected line broadening due to electron collisions, 
C 4 = 0. The results are listed in Tables Q, [| |J and @. Most 
of the lines are 'fictitious' lines ('even' wavelength entries in 
the tables), allowing a systematic study of the granulation ef- 
fect as a function of the line parameters. The tables also con- 
tain a number of 'real' lines (table entries with 'uneven' wave- 
length) for Li, O, and Si, which have alrea dy bee n used with de- 
tailed abundance studies by Wedemeyer (2001) and Holweger 
(2001). In the latter work, the granulation abundance correc- 
tions were obtained by interpolation in the list of fictitious lines 
(Tables [j],p]), which covers a sufficient range of the relevant line 
parameters. 

5. 1. Atoms with high ionization potential and ions 



Atoms with high ionization potential like C I, N I, and O I, as 
well as singly ionized elements are representative of species in 
the main ionization stage. Neutral nitrogen (ionization potential 
X = 14.5 eV) can be taken as the prototype of these species. 
Unlike C and O, N is not significantly affected by molecule 
formation or changes in the degree of ionization. 

The corrections A gran as a function of excitation energy, 
Ei, are displayed in Fig. |]for our tabulated sample of weak 
fictitious lines (A 5500 A, W s» 5 mA) of these species. 
Obviously, the corrections for high-excitation lines of the 'ma- 
jority species' are similar for all elements, and can amount to 
-0.1 dex for weak lines with Ei *S 10 eV. This conclusion is 
independent of the assumed velocity field which is irrelevant 
here. A mn _, 



A gran b for lines as weak as 5 mA (see ta- 



bles). 
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Fig. 5. Logarithmic granulation abundance corrections A gran 
as a function of excitation energy Ei for C I, N I, O I, and a 
number of singly ionized elements, representative of species 
in the main ionization stage. The plot is based on a subset of 
weak fictitious lines with A = 5 500 A, W\ ~ 5 mA, taken 
from Tables jjj || @. The upper symbol at E t = 12 eV for N I 
refers to A 10 000 A. 
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Fig. 6. Granulation abundance corrections, A gran , for the first 
8 fictitious nitrogen lines listed in Table [j], compared with test 
results obtained by increasing the line strength by factors of 5 
and 15, respectively. 



5.1 .1 . Dependence on line strength 

The dependence of the granulation abundance corrections on 
line strength was investigated by test calculations for two ad- 
ditional sets of stronger N I lines. The two set of lines were 
obtained from the set of N I lines (A 5 500 A) listed in Table |l] 
by increasing their equivalent widths, Wm, by a factor of 5 
and 15, respectively. Thus, the test lines in first set of have 
Win ~ 25 mA, those in the second set Win ~ 80 mA. 
According to the results shown in Fig. [| the granulation abun- 
dance corrections A gran are systematically more positive for 
the stronger, otherwise identical lines. For Ei = 12 eV the dif- 
ference in A gran between the lines with Wm ~ 80 mA and 
5 mA amounts to +0.07 dex. 




Fig. 7. Granulation abundance corrections, A gran , for the first 
8 fictitious nitrogen lines listed in Table [lj compared with test 
results obtained by moving the weak lines from A = 5 500 A to 
different spectral regions, centered at A = 4 000, 8 000, 10 000, 
12 000, and 16 000 A, respectively, keeping the reduced equiv- 
alent width W/X fixed (Wm ~ A/5 500 • 5 mA). 



5.1 .2. Dependence on wavelength 

In order to get an idea about the dependence of the granula- 
tion abundance corrections on the wavelength of the spectral 
lines, we also performed test calculations for five additional 
sets of N I lines, obtained by shifting the weak N I lines listed 
in Table @ from A = 5 500 A to A = 4 000, 8 000, 10 000, 
12 000, and 16 000 A, respectively, enforcing the same reduced 
equivalent width, W/X, i.e. Wm ~ A/5 500 ■ 5 mA. 

The results shown in Fig. ^ indicate that the granulation 
abundance corrections A gran indeed depend on wavelength, in 
the sense that lines in the blue part of the spectrum show sys- 
tematically more negative corrections than those in the red part. 
This general trend is evident over the whole wavelength range, 
apart from a slight reversal between A 8 000 and 10 000 A 
which is related to the Paschen jump at A 8208 A. For Ei = 
12 eV the difference in A gran between a weak line at A 4 000 
and 10 000 A amounts to -0.1 dex. 



5.2. Atoms with low ionization potential 

Elements with low ionization potential are predominantly 
singly ionized in the solar atmosphere. The particle density of 
the remaining neutral atoms of such elements is particularly 
sensitive to temperature variations. Molecules fall into the same 
category. Neutral iron (ionization potential \ = 7.9 eV) can 
be taken as representative for these species. According to our 
models, the fraction of neutral iron is nowhere exceeding 26%; 
it shows a maximum near logrssoo = —3 and has dropped to 
less than 3.5% at logr 55 oo = 0. 

The corrections A gran for a sample of weak fictitious lines 
(A 5500 A, W ~ 5 mA) of these 'minority species' are dis- 
played in Fig. H as a function of excitation energy, E^ We note 
that the largest granulation effects are found for low-excitation 
lines of neutral atoms with an ionization potential between 6 
and 8 eV. The magnitude of the correction depends on the el- 
ement and is as large as -0.3 dex for Til, Ej = eV. Again, 
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Fig. 8. Logarithmic granulation abundance corrections A gran 
as a function of excitation energy Ei for neutral atoms of pre- 
dominantly ionized elements, and for molecules CN and MgH. 
The plot is based on a subset of weak fictitious lines with 
A = 5 500 A, W x « 5 mA, taken from Tables § (§ § f§ 



Wavelength dependence of A aran B , WA fixed 
0.05 P ' ' ' 1 ' ' ' 




-0.30 E 
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Fig. 10. Granulation abundance corrections, A gran , for the 7 
fictitious neutral iron lines listed in Table ||, compared with test 
results obtained by moving the weak lines from A = 5 500 A to 
different spectral regions, centered at A = 4 000, 8 000, 10 000, 
12 000, and 16 000 A, respectively, keeping the reduced equiv- 
alent width W/X fixed (W 1D w A/5 500 • 5 mA). 




Ei [eV] 



Fig. 9. Granulation abundance corrections, A gran , for the 7 
fictitious neutral iron lines listed in Table || compared with test 
results obtained by increasing the line strength by factors of 5 
and 10, respectively. 



we point out that this result does not depend on the assumed 
non-thermal velocity field. 



It is interesting to note that Asplund ( 2000a I, using 3D 
models, finds corrections in the opposite direction for the 'real' 
Si I lines listed in Table [| We emphasize that this is not an 
indication of fundamentally different "granulation corrections" 
in 2D and 3D. The corrections given by Asplund with respect 



to the Holweger-Miiller (1974) model include not only the ef- 



fect of temperature inhomogeneities but also the influence of 
the different mean stratifications. Obviously, the difference in 
the mean structure is the more important of the two opposing 
effects in this case. Indeed, we find from our own models a to- 
tal correction of -0.05 dex with respect to the Holweger-Miiller 
model, compared to a the tabulated "granulation correction" of 
+0.02 dex, explaining the apparently inconsistent results. 



5.2.1 . Dependence on line strength 

As for N I, the dependence of the granulation abundance cor- 
rections on line strength was investigated by test calculations 
for two additional sets of stronger Fe I lines, obtained from the 
set of Fe I lines listed in Table |] by increasing their equiva- 
lent widths, IVid, by a factor of 5 and 10, respectively. The 
test lines in first set of have then Wiu ~ 35 mA, those in the 
second set W±d ~ 70 mA. The results are shown in Fig. ^[ 
Qualitatively, the granulation abundance corrections A gran for 
Fe I show the same dependence on line strength as already seen 
for N I: they are systematically more positive for the stronger, 
otherwise identical lines. 

The corrections for Ei = and 1 eV, Wir, ~ 70 mA, 
shown in Fig. ^|are certainly too small (too positive), because a 
significant part of the line absorption comes from layers lo- 
cated outside the model atmosphere on which the line for- 
mation calculations are based (T5500 < 10~ 5 ). This situation 
results in a stronger underestimation of the equivalent width 
for the inhomogeneous case. Hence, A gran should actually 
be more negative, and we can expect abundance corrections 
A gra n < — 0.1 dex even for the stronger ground state lines. 

5.2.2. Dependence on wavelength 

We have performed test calculations for five additional sets of 
Fe I lines, obtained by shifting the standard Fe I lines listed in 
Table | from A = 5 500 A to A = 4 000, 8 000, 10 000, 12 000, 
and 16 000 A, respectively, enforcing the same reduced equiv- 
alent width, W/X, i.e. W 1D w A/5 500 • 5 mA. 

The results shown in Fig. [jjo] indicate that the wavelength 
dependence of the granulation abundance corrections A gran for 
Fe I has the opposite sign as in the case of N I. Lines in the blue 
part of the spectrum tend to show smaller (less negative) cor- 
rections than those at longer wavelengths. This general trend is 
clearly evident from A 4000 to 8000 A; at longer wavelengths 
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the variation of A gran with A is much reduced and no longer 
strictly mono tonic (for an explanation see Paper II). 

We have to point out that the results for A 16 000 A are un- 
certain, because a significant part of the line absorption comes 
from layers above T16000 < 10 -5 , so the line formation region 
is not entirely covered by our model atmosphere. 



6. Conclusions 

In an effort to estimate the - hitherto largely unknown - ef- 
fects of photospheric temperature fluctuations on spectroscopic 
abundance determinations, we have carried out numerical sim- 
ulations of the problem of LTE line formation in the inhomo- 
geneous solar photosphere based on detailed 2-dimensional ra- 
diation hydrodynamics models of the convective surface layers 
of the Sun. 

For a variety of spectral lines of different elements, we have 
computed synthetic line profiles from the inhomogeneous 2D 
hydrodynamical atmosphere and from the corresponding ID 
plane-parallel model, respectively. By means of a strictly dif- 
ferential 1D/2D comparison of the emergent equivalent widths, 
we have derived so called "granulation abundance corrections" 
for the individual lines, which have to be applied to standard 
abundance determinations based on homogeneous ID model 
atmospheres in order to correct for the influence of the pho- 
tospheric temperature fluctuations. The 'classical' problem of 
finding the most appropriate mean vertical temperature stratifi- 
cation is not addressed here. 

Using the concept of "fictitious" spectral lines, we were 
able to investigate systematically the dependence of the "granu- 
lation abundance corrections" on the basic line parameters like 
excitation potential, line strength, wavelength, element and ion- 
ization stage. For many lines of practical relevance, it should be 
possible to estimate the magnitude of the abundance correction 
by interpolation in the graphs and tables provided in this pa- 
per. This approach may often be an acceptable alternative to 
a detailed fitting of individual line profiles based on hydrody- 
namical simulations. 

In general, we find a line strengthening in the presence of 
temperature inhomogeneities, implying mostly negative "gran- 
ulation abundance corrections", i.e. standard analysis based on 
plane-parallel atmospheres tends to overestimate abundances. 
The physical reason for the line strengthening is primarily the 
non-linear temperature dependence of the line opacity due to 
thermal ionization and excitation, as demonstrated in Paper II. 

One remarkable result is that all lines of our sample with an 
excitation potential around Ei w 5 eV are practically insensi- 
tive to granulation effects, regardless of element and ionization 
stage. Moderate granulation corrections (A gran w —0.1 dex) 
are found for weak, high-excitation lines (Ei 10 eV) of ions 
and atoms with high ionization potential like N I. The largest 
corrections are found for ground state lines (Ei = eV) of 
neutral atoms with an ionization potential between 6 and 8 eV 
like Mgl, Cal, Til, Fel, amounting to A glan sw —0.3 dex in 
the case of Ti I. 

For given excitation potential, the granulation corrections 
are systematically more positive for stronger lines. The wave- 



length dependence of A gran , however, depends on the type of 
line ('majority' or 'minority' species). 

We recall that the corrections A gran account only for the 
influence of the temperature and pressure fluctuations relative 
to a given mean structure. We have suppressed possible effects 
related to 1D/2D differences in the details of the non-thermal 
velocity field. Such additional effects, which apply to partly 
saturated lines only, should be small, however, since they are 
largely eliminated by the proper choice of the microturbulence 
parameter £ m i c . 

We point out that the magnitude of the "granulation abun- 
dance corrections" derived in this work should be considered 
as upper limits (\A gl . an \{2D, LTE) > |A gran |(3£>, NLTE)), 
because (i) possible NLTE-effects, which were ignored in this 
study, tend to reduce the fluctuations of the line opacity (e.g. 
Kiselman |1997|, Cayrel & Steffen fZOM Asplund |2000b[), and 



(ii) the amplitude of the temperature fluctuations is systemat- 
ically overestimated in 2D relative to 3D convection models 
(see Fig. Q). The large negative corrections found for the low- 
excitation lines of atoms with low ionization potential are the 
result of the strongly increasing amplitude of the temperature 
fluctuations towards the upper photosphere as predicted by our 
present 2D hydrodynamical simulations. In a more realistic 3D 
atmosphere, this feature will therefore be less pronounced. 

Even though the results are still somewhat uncertain, the 
present study should be helpful in providing upper bounds for 
possible errors of spectroscopic abundance analyses, and for 
selecting spectral lines which are as insensitive as possible to 
the effects of photospheric temperature inhomogeneities. 

For a physical explanation of the numerical results obtained 
in this work, a detailed investigation of the process of LTE line 
formation in an inhomogeneous stellar atmosphere based on 
the analysis of the transfer equation and the evaluation of line 
depression contribution functions is presented in the second pa- 
per of this series (Paper II). With the help of this formalism, we 
can actually tell how the different atmospheric layers contribute 
to the granulation correction for any particular spectral line. 

In the near future, we plan to extend the computation of 
abundance corrections to full 3D simulations for the Sun and 
other types of stars where convection extends into higher pho- 
tospheric layers and the spectroscopic granulation effects may 
well be even larger. 
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(5. /O/) 


-0.078 
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Oi 
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0.0 
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5.076 


(5.069) 


+0.015 


(+0.016) 


0.003 


Oi 


5500.0 


2.0 


5.088 


4.847 


(4.838) 


+0.022 


(+0.023) 


0.007 


Oi 
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4.0 


5.053 
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+0.011 
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Oi 
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Oi 
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(5.725) 


-0.045 


(-0.042) 


0.009 


1 
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1 
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(+0.011) 
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Oi 
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5.408 


6.375 


(6.336) 


-0.081 


(-0.078) 


0.014 


Oi 


6300.3 


0.00 


4.479 


4.329 


(4.324) 


+0.015 


(+0.016) 


0.003 


Oi 


7771.9 


9.15 
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81.184 


(79.130) 


+0.013 


(+0.039) 


0.015 


Oi 
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9.15 
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(64.988) 


+0.007 


(+0.031) 


0.013 


Oi 


7775.4 


9.15 


51.250 


^1 111 
Jl.Z/j 




-0.000 


(+0.021) 


0.012 


Oi 


9265.9 


10.74 


34.068 


35.492 


(35.062) 


-0.026 


(-0.018) 


0.009 


Oi 


11302.4 


10.74 


14.074 


14.947 


(14.807) 


-0.033 


(-0.028) 


0.008 


Oi 


13164.9 


10.99 


17.303 


17.894 


(17.705) 


-0.017 


(-0.013) 


0.006 



Table 1. Granulation abundance corrections for Li I, C I, N I, and O I, derived from a sample of representative 'fictitious' lines 
(even wavelengths) and some selected 'real' lines commonly used in the analysis of late type stars. E\ is the excitation potential 
of the line's lower energy level, W2D is the mean equivalent width computed from 1 1 selected 2D snapshots, replacing the 
hydrodynamical velocity field by a depth-independent mircroturbulence of £ m i c = 1 km/s; Wid is the mean equivalent width 
computed from the corresponding ID models (£ m i c = 1 km/s). A gran is the resulting granulation abundance correction derived 
from the ratio /Win, including curve-of-growth effects. For reference, we give in parentheses the results obtained when the 
2D line profiles are computed with the hydrodynamical velocity field of the respective snapshot; these corrections, however, are 
not recommended (see text for details), a is the standard deviation of the abundance corrections derived from the 1 1 individual 
snapshots, measuring the statistical uncertainty of the tabulated A gran . 
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0.012 
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5.709 
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(5.439) 


+0.013 
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Si I 
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38.063 


36.710 


(36.131) 


+0.022 


(+0.031) 
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Si I 
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90.285 


88.150 


(85.562) 


+0.023 


(+0.051) 


0.006 


Si I 
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56.938 


55.137 
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(4.553) 


+0.025 
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4.921 
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(4.816) 


+0.008 


(+0.010) 
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5.057 


5.234 
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-0.016 


(-0.012) 
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(5.415) 


-0.031 


(-0.026) 
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fx 
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(-0.074) 


0.014 


Sill 
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11.0 


5.506 


6.633 


(6.476) 


-0.104 


(-0.091) 


0.017 


Sin 


6347.1 


8.12 


47.386 


48.453 


(46.572) 


-0.021 


(+0.016) 


0.013 


Sin 


6371.4 


8.12 


33.817 


34.876 


(33.622) 


-0.025 


(+0.005) 


0.011 



Table 2. Same as Table [j], but for a sample of Na I, Mg I, Mg II, Si I, and Si II lines. 
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-0.005 


z / \ /www 

(-0.000) 
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Fel 


5500.0 


3.0 
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0.0 
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(4.813) 


+0.028 


(+0.031) 


0.007 


Fell 
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2.0 


5.076 
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(4.761) 


+0.026 


(+0.029) 


0.007 


Fell 


5500.0 


3.0 


5.040 


4.821 


(4.777) 


+0.020 


(+0.025) 


0.007 


Fell 


5500.0 


4.0 


5.031 


4.891 


(4.838) 


+0.013 


(+0.018) 


0.007 


Fell 


5500.0 


5.0 


5.058 
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Fe II 


5500.0 
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5.224 


5.618 


(5.475) 


-0.037 


(-0.024) 


0.009 


Fell 
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3.0 


5.040 


4.821 


(4.777) 


+0.020 
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88.180 
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+0.028 


(+0.032) 


0.008 
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3.0 


61.639 
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+0.030 


(+0.139) 


0.013 


Fell 


5500.0 


3.0 


60.639 


58.560 


(53.935) 


+0.041 


(+0.134) 


0.014 


Fell 


6500.0 


3.0 


59.936 


57.334 


(52.999) 


+0.048 


(+0.130) 


0.014 


Fell 


7500.0 


3.0 


59.520 


56.536 


(52.458) 


+0.053 


(+0.126) 


0.015 



Table 3. Same as Table [j], but for a sample of fictitious S I, Ca I, Ca II, Ti I, Ti II, Fe I, and Fe II lines. 
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0.061 



Table 4. Same as Table [j], but for a sample of fictitious Sr I, Sr II, Ba II, CN, and MgH lines. 



